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The fusion excitation functions have been measured with rather good accuracy for ^^S+'^Zr and ^^S+'^Zr 
near and below the Coulomb barrier. The sub-barrier cross sections for ''^S+'^Zr are much larger compared 
with ^"S+'^Zr. Semi-classical coupled-channels calculations including two-phonon excitations are capable to 
describe sub-barrier enhancement only for '"S+'^Zr. The remaining disagreement for ^^S-l-'^Zr comes from the 
positive 2-value intermediate neutron transfers in this system. The comparison with '"'Ca+'*Zr suggests that 
couplings to the positive 2- value neutron transfer channels may play a role in the sub-barrier fusion enhance- 
ment. A rather simple model calculation taking neutron transfers into account is proposed to overcome the 
discrepancies of ''^S+'^Zr. 

PACS numbers: 25.70.Gh, 25.70.Jj, 24.10.Eq 



I. INTRODUCTION 

The heavy-ion fusion reactions in the low-energy range 
near and below the Coulomb barrier have been the subject 
of extensive experimental and theoretical efforts in the past 
decades [l,-t3J. Beside the fact that the questions of the pos- 
sible occurence of unexpected phenomena, such as breakup 
effects on the fusion reactions at near barrier energies 13|], are 
still unresolved, one has still to understand better the role of 
neutron transfers in the fusion process Jj, |51]. For instance, 
effects of neutron-rich projectiles on the formation of super- 
heavy elements (SHE) |6], especially with the development of 
newly available Radioactive Ion Beams (RIB) facilities need 
to be clarified as well as fusion hindrance at extremely low 
energies that remain among the most interesting open ques- 
tions in the nuclear astrophysics domain |7]. Fusion enhance- 
ment below the Coulomb barrier is one of the most studied 
phenomena and, measurements of fusion barrier distributions 
have been widely performed to investigate the mutual impor- 
tance of both the nuclear structure and dynamical process ef- 
fects on the sub-barrier fusion enhancement Il8l- [l4ll . 

Coupled-channels (CC) calculations have been used to de- 
scribe the reactions in this energy range theoretically (see for 
example Refs. 1 1, 8] and references therein). Fusion enhance- 
ment due to the static deformations and surface vibrations of 
the nuclei has been well described in the coupled-channels 
calculations lH-H. 

The influence of the neutron transfer channels on sub- 
baiTier fusion process 111-171 is not yet fully understood. 
During the last 20 years a large number of experimental 
and theoretical investigations were undertaken to study the 
neutron-transfer mechanisms in competition with the fusion 
process. Stelson et al. iflSi - UTIl proposed an original scenario 
that uses an empirical method involving a sequential transfer 



of several neutrons between the reactants. This multineutron 
transfer process is capable to initiate fusion at large internu- 
clear distances and will smooth the fusion barrier distributions 
(with larger width) with lower energy thresholds. This "shift" 
effect coiTesponds to the energy window for which the nu- 
clei are allowed to come sufficiently close together for neu- 
trons to flow freely between the target and projectile. As 
a consequence, this will reduce the effective barrier and en- 
hance the fusion cross sections at sub-barrier energies. Fol- 
lowing this idea, Rowley et al. [18j have first used a simple 
phenomenological model that simulates coupling to neutron 
transfer channels with a parametrized coupling matrix. Later 
on, Zagrebaev [19] proposed another semiclassical theoreti- 
cal model that has been successfully used to reproduce the 
sub-barrier fusion enhancement of the old '*"Ca-H'^''Zr reaction 
data of Ref. [11] by including the intermediate positive Q- 
value neutron transfer channels in the CC calculations. 

The failure of the CC calculations including only the 
couplings to the inelastic excitations indicates that couplings 
to neutron transfer channels may play a key role in the fusion 
dynamics near the barrier for medium-heavy systems such 
as ^OCa+^O-^^Zr d M, ^^S+'^O'^^Zr M, and ^ONe+^O-'^^Zr 



|22|1, the last two reactions being studied by measurements of 
large-angle quasielastic scattering. Previous measurement of 
quasielastic scatterings of ■'^S-i-'^'''''''Zr were also undertaken 
at backward angles near the barrier Il23;l : their analysis gave 
indication that positive Q-value neutron transfer channels 
should be included in the coupling scheme. Up to now 
no fusion data exist neither for ^"Ne-i-^"'^^Zr nor for the 
32(g^_90,962j- reactions, it will be interesting to measure their 
fusion excitation functions. In order to disentangle the 
possible effect of positive Q-value neutron transfer couplings 
we decided to investigate the two last systems. We report 
here about the measurement of near- and sub-barrier fusion 



excitation functions of ■'^S+^^Zr and ^^S+^*'Zr performed 
with small energy steps and good statistics accuracy. 

Our research will focus on the role of neutron transfers be- 
tween the colliding nuclei as a mechanism that enhances the 
fusion cross sections at sub-barrier energies. This paper is or- 
ganized as follows: Sec. II presents the experimental setup 
and details on the measurements. Results of the analysis of 
the experimental data are given in Sec. III. Their discussion 
is finally proposed in Sec. IV in the framework of compar- 
isons with semiclassical coupled-channels calculations before 
a short summary of the Sec. V. 



II. EXPERIMENTAL PROCEDURES 

The experiment was performed at the HI- 13 tandem accel- 
erator of CIAE, Beijing. The collimated ^^S (Q - 10^ charge 
state) beam was used to bombard the zirconium oxide targets. 
The beam intensity was stabilized in the 2-20 pnA range in 
order to minimize the pile-up for each of the bombarding en- 
ergies. The 3 mm diameter (98.87% enriched) "^"ZrOa and 
(86.4% enriched) ^^ZrOa targets were both 50 fig/cm^ thick 
and evaporated onto 15 (xg/cm^ carbon foil backings. The 
beam energies varied over the range Ei-^h = 100-130 MeV for 
^^S+'^°Zt and £iab = 95-130 MeV for ^^S-n'^^Zr (in 1.33 Mev 
steps at the higher and 0.67 MeV at the lower energies) and 
changed only downwards starting at fiiab = 130 MeV in order 
to reduce the magnetic hysteresis for both targets. The target 
chamber contains four silicon detectors at = +25 ° symmet- 
rically (right/left and up/down) with respect to the beam axis 
in the forward direction in order to monitor the beam optics 
(Rutherford scattering) and to provide an absolute normaliza- 
tion of the fusion cross sections. 

The fused evaporation residues (ER) concentrated to within 
a few degrees of the incident beam direction were separated 
from the incident beam (see Fig. 1) by an electrostatic deflec- 
tor which design is pretty much similar to the experimental 
setup in Legnaro |24]. It consists of two pairs of electrodes 
followed by an E-TOF arrangement with of a microchannel 
plate (MCP) detector coupled to a Si(Au) surface barrier de- 
tector Two-dimensional plots of the data were used to cleanly 
separate the ER's from the beam-like products (BLP). A typ- 
ical example of the time-of-flight versus energy spectrum for 
^^S-H'^^'Zr measured at Ei^b = 130 MeV and = 2 ° is shown 
in Fig. 1. The electrostatic deflector could be rotated about 
the target position in the horizontal plane to measure the ER 
angular distributions. 

The particles coming from the target were selected before 
entering the fields by an entrance collimator of 3 mm diam- 
eter, corresponding to a A0 = + 0.57 ° opening. A lOjUg/cm^ 
thick carbon foil clung to the collimator was used to reset the 
atomic charge state distribution on the ion path. The colli- 
mator of the MCP defined the solid angle of the electrostatic 
deflector as being approximately AQ. = 0.3 msr 

ER angular distributions were measured in the range = 
-4 ° to = 10 ° with step A0 = 1 ° at three beam energies (E\.^b 
- 100, 115, and 130 MeV) for both systems. The angular dis- 
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FIG. 1: (Color online) Two-dimensional plot E-TOF of the events, 
following beam separation, taken at iii^b =130 MeV and at 2° for 
the '-S+'^Zr reaction. Two groups of particles [beam-like parti- 
cles(BLP) and evaporation residues (ER)] are indicated. 



tributions were found to be symmetrical about = °, as ex- 
pected. Their typical shape did not change appreciably with 
the beam energy. These combined angular distributions and 
the double Gaussian fits were used to obtain the fusion cross 
sections. At each energy the number of ER events was nor- 
malized to the Rutherford scattering rates counted by the mon- 
itor detectors. For the most of energy points, only differential 
cross sections were measured at 9 -2°, from the obtained val- 
ues, the total ER cross sections were deduced. Using the solid 
angles, the = 2 °-to-total ratios and the measured transmis- 
sion efficiencies, these ER yields were transformed into total 
cross sections. Since fission of the compound nucleus can be 
neglected for both systems, the measured cross sections were 
taken as complete fusion cross section Of. 

The transmission efficiencies and the relevant voltages used 
to deflect the ER were calibrated by the '^^Ba beam scattered 
by the ^"Zr target at small angles and at the corresponding 
energies with the fusion evaporation residues. It was found 
that the defocusing effect of the deflection voltage reduces 
the transmission from unity to 0.60±0.06. Additional system- 
atic errors come from the geometrical solid angle uncertain- 
ties, the angular distribution integrations, and the transmission 
measurements. Altogether these contributions sum up to a + 
15% value for systematic errors. 



III. EXPERIMENTAL RESULTS 

The measured ER excitation functions for the two systems 
are shown in Fig. 2, where the energy scale is corrected for 
the target thickness. The statistical errors shown in the figure 
do not exceed the symbol size for most of the experimental 
points. They are + 0.8% for both the high-energy and the 
intermediate-energy points and increase to + 23% for the low- 
energy points. The ER cross sections are listed in Tables I and 
II for both reactions. 



TABLE I; Experimental ER cross sections for ^"S+'^Zr. 



£,„. (MeV) 


o-ER (mb) £, 


m. (MeV) 


o-ER (mb) 


95.3 


623.51+4.00 


82.7 


142.10+1.37 


94.2 


582.98+4.63 


82.2 


122.84+1.08 


93.3 


528.80+3.93 


81.7 


105.94+1.10 


92.2 


535.48+4.70 


81.2 


91.59+0.92 


91.2 


430.00+3.37 


80.7 


77.96+0.80 


90.7 


442.29+3.03 


80.2 


63.69+0.62 


90.2 


429.44+3.58 


79.7 


54.96+0.55 


89.7 


363.05+3.22 


79.2 


40.22+0.39 


89.2 


357.17+2.84 


78.7 


32.89+0.34 


88.7 


342.70+2.94 


78.2 


23.41+0.35 


88.2 


322.27+3.07 


77.7 


17.46+0.35 


87.7 


309.41+2.80 


77.2 


10.58+0.42 


87.2 


298.07+2.85 


76.7 


5.89+0.41 


86.7 


277.99+2.29 


76.2 


3.66+0.37 


86.2 


246.27+2.81 


75.7 


2.07+0.25 


85.7 


244.09+2.02 


75.2 


1.09+0.16 


85.2 


224.93+1.95 


74.7 


0.57+0.10 


84.7 


190.06+1.64 


74.2 


0.31+0.06 


84.2 


189.70+1.75 


73.7 


0.18+0.04 


83.7 


167.85+1.25 


73.2 


0.12+0.03 


83.2 


143.85+1.38 







TABLE IL Experimental ER cross sections for ""^S+'^Zr. 



£c.m. (MeV) 


o-ER (mb) £, 


.m. (MeV) 


o-ER (mb) 


97.1 


683.96+3.93 


82.1 


161.35+1.34 


95.9 


644.56+7.76 


81.8 


145.06+1.51 


95.0 


601.59+3.70 


81.1 


130.62+1.09 


93.9 


591.22+4.33 


80.8 


117.57+1.21 


93.0 


505.50+3.00 


80.1 


99.28+0.85 


92.3 


519.80+4.10 


79.7 


90.53+0.83 


92.0 


510.25+3.57 


79.1 


72.13+0.66 


91.3 


439.19+3.38 


78.7 


66.85+0.70 


90.9 


448.78+3.14 


78.1 


50.41+0.43 


90.3 


381.53+2.95 


77.7 


46.28+0.51 


89.9 


424.93+3.13 


77.0 


37.01+0.31 


89.3 


382.96+3.02 


76.7 


31.81+0.34 


88.9 


393.91+2.55 


76.0 


23.82+0.03 


88.3 


354.74+2.11 


75.7 


20.18+0.03 


88.0 


338.00+2.74 


75.0 


15.32+0.23 


87.2 


314.40+2.44 


74.6 


13.12+0.26 


86.9 


313.05+2.37 


74.0 


8.95+0.36 


86.2 


264.38+2.37 


73.6 


5.53+0.39 


85.8 


282.11+2.51 


73.0 


2.95+0.30 


85.2 


244.76+1.78 


72.5 


1.94+0.23 


84.8 


231.95+2.18 


71.9 


0.79+0.12 


84.2 


222.32+1.83 


71.4 


0.69+0.12 


83.8 


231.63+1.93 


70.9 


0.37+0.07 


83.2 


182.87+1.85 


70.4 


0.23+0.05 


82.8 


179.30+1.82 


69.9 


0.09+0.02 
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FIG. 2: Experimental ER excitation functions of ''^S+"'Zr (hollow 
circles) and ^^S+'*Zr (solid circles) as a function of the center-of- 
mass energy. The error bars represent purely statistical uncertainties. 
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FIG. 3: Plot of the ER excitation functions on a reduced-energy scale 
for the two systems ^^^S+'*Zr (filled circles) measured in the present 
work and '^''Ca+'^Zr (open circles). The data for '"'Ca+'*Zr are taken 
from Ref. fill]. 



The comparison of the ER excitation functions of ^^^S+''*Zr 
(present work) and '^^Ca+^^Zr 1 1 1] is made easier when cross 
sections are plotted in a reduced-energy scale as shown in Fig. 
3. One observes that the two systems display very similar be- 
haviors on the whole energy range despite relatively large dis- 
crepancies. These discrepancies are mainly due to larger un- 
certainties in the present data arising from larger backgrounds 
in the spectra. It is interesting to notice that both system reac- 
tants have very similar nuclear structures as well as neutron- 
transfer properties. This behavior, already discussed in our 
previous investigation of ^^S-i-^'^-^^Zr quasielastic barrier dis- 
tributions ['23'], indicates that the positive Q-value neutron 
transfers strongly enhance the fusion cross sections at sub- 
barrier energies. This experimental observation will be con- 
firmed by the semiclassical coupled-channels calculations as 
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FIG. 4: ER excitation function of ^^^S+'^Zr. The open circles are 
the experimental data. The dotted, dashed, and solid lines represent 
the CCFULL calculations without coupling and with the one-phonon 
and two-phonon couplings, respectively (see text for details). The 
arrow indicates the position of the Coulomb barrier for ^^^S+'"'Zr. 



FIG. 5: ER excitation functions of ''^S+'^Zr. The solid circles are 
the experimental data. The dotted, dashed and dash-point lines are 
the CCFULL calculations with no coupling, one-phonon and two- 
phonon couplings, respectively. The solid line is the calculation tak- 
ing into account the neutron transfers by NTFus code. The arrow 
indicates the position of the Coulomb barrier for '-S+'^Zr. 



discussed in the following Section. 



IV. DISCUSSION: SEMICLASSICAL 
COUPLED-CHANNELS CALCULATIONS 

The ER excitation functions of ^^S+'^"Zr and ^^S+^^7x have 
been compared with the results of CC calculations performed 
with the CCFULL code ll26ll using the Akyiiz-Winther nuclear 
potential 1271]. 

The relevant informations on the low-lying excitations of 
^^S, '"Zr, and ^*Zr can be seen in Table IIL The quadrupole 
vibrations of both ^"Zr and ^^Zr nuclei are weak in energy: 
in fact, they lie at comparable energies. With this potential 
the CCFULL barriers were found to be at Vb = 81.2 MeV for 
^^S+^°Zr and at V^ = 80.1 MeV for ^^S+'^^Zx, respectively. 
These values are fully consistent with what was found for the 
quasielastic barriers previously measured by our group ll23ll . 
In a first approximation, these values can be considered as av- 
erage values between the barrier heights for the nose-to-nose 
configuration and for the side-to-side configuration l25[l . 

Figures 4 and 5 show the comparisons of the experimental 
ER excitation functions and the CCFULL calculations with 
one- and two-phonon couplings and without coupling for the 
32(g_|_90,%2j- fusion reactions. Semiclassical coupled-channels 
calculations have also been performed by use of the new code 
NTFus 1 28] where the proximity potential is adopted. The 
new oriented object code has been constructed in the frame- 
work of the Zagrebaev model ] 25 ] and implemented in C^^ by 
using the compiler of ROOT ]29]. More technical details and 
related discussions about the NTFus code will be illustrated 
elsewhere ]28]. 

First we propose to present the calculations for "'^S-n'^'^Zr 
that are displayed in Fig. 4. They were performed without 
taking the neutron transfers into account. The two-phonon 
coupling CCFULL calculations are quite satisfactory. The 



calculation (solid line) reproduces the data below and above 
the barrier Vb (arrow in Fig. 4). On the other hand, the cal- 
culations fail for ^^S-i-''*'Zr with large discrepancies occuring 
mainly at energies below the barrier V^ (arrow in Fig. 5). 
One- and two-phonon (both shown in Fig. 5) excitations 
in ^''Zr, of both quadrupole and octupole natures, bring ad- 
ditional but not sufficient enhancements. Finally, we tried 
also "three-phonon" couplings, but no further improvement 
could be reached. Anyway, additional couplings in ^^S-i-'^*Zr, 
which might give rise to lower-energy barriers, are simply 
not present in the coupling scheme. Similar conclusions have 
been obtained for the "^"Ca-H^'^Zr reaction Ill4ll . 

When we choose to take into account the neutron transfers, 
the fusion excitation function can be derived using the follow- 
ing formula 1191 : 



r 1 v-i r^°**' 

Ti(E) = /(B) ttX a,(E,l,Q) 

J ^tr ■^ J-E 

xPhw(B, E + Q, OdQdB, 



and 



nJf 



la- 



(Tfus(E)^—-yQl+\)Ti{E), 



(1) 



(2) 



where Ti(E) is the transmission, E the energy at the center-of- 
mass, f{E) the normalized barrier distribution function, I the 
momentum and l„ the critical momentum calculated where 
there is no coupling (well above the barrier). a\SE,l,Q) and 
2o(k) are the probability and the g-value for the transfer of 
k neutrons, 1/N„ is the normalization of the total probability 
taking into account the neutron transfers. 

The calculation with the neutron transfer effect is per- 
formed up to the channel -i-4n (k=4). No more visible effect 



TABLE III: Excitation energies Ex, spin and parities A", and defor- 
mation parameters/?,! for ^^^S and '"-'^Zr. 



Nucleus 


£x (MeV) 


r Ai 


3^8 


2.230 


2+ 0.32 




5.006 


3- 0.40 


"OZr 


2.186 


2+ 0.09 




2.748 


3- 0.22 


"'Zr 


1.751 


2+ 0.08 




1.897 


3- 0.27 



TABLE IV: 2-value in MeV for neutron pickup transfer channels 
from ground state to ground state for the '-S+'^'^Zr systems. 



System +ln +2n +3n +4n 
32s+*Zr -3.33 -1.229 -6.59) -6.319 
^^S-h'^Zr 0.788 5.737 4.508 7.655 



can be obtained by using -i-5n and +6n channels. More de- 
tails of the calculation procedures and of the description of 
the NTFus code [28], itself, will be given in a forthcoming 
Brief Report. The Q-values for the calculation (solid line in 
Fig. 5) are given in Table IV. As we can see in Fig. 5, the 
dash-point line (without the neutron transfers) does not at all 
describe the data at the sub-barrier energies. In contrast, the 
solid line taking into account the neutron transfers is able to fit 
the data reasonably well. As expected, the correction applied 
on the calculation at sub-barrier energies by the Zagrebaev 
model III9II25I1 enhances the cross sections further. Moreover, 
it allows a fairly good description of the present experimen- 
tal data showing the strong effect of neutron transfers for the 
sub-barrier fusion of ■'^S-i-^^Zr 

Figures 6 and 7 show the experimental barrier distributions 
from fusion and quasi-elastic scattering and the corresponding 
CCFULL calculations for the two systems. The fusion barrier 
distributions for the two systems have been obtained by dou- 
ble differentiation Ecrfu~i vs energy using the three-point differ- 
ence formula 1 1 IL The quasi-elastic barrier distributions are 
taken from Ref . ||23I1 . It is very interesting to note that for both 
reactions the experimental quasi-elastic barrier distributions 
and the experimental fusion barrier distributions are strikingly 
similar The large fluctuations occurs in the barrier distribu- 
tions for Ec.m. >80 MeV is not enough due to the measuring 
accuracy of the ER cross sections. For ^^S-i-^^Zr, the over- 
all trends of the experimental barrier distributions are roughly 
consistent with the CCFULL calculation considering the two- 
phonon coupling. While for ^-SH-''^Zr, the experimental bar- 
rier distributions are wide and show a low-energy tail extend- 
ing to the lowest energies compared with ^^S+''"Zr and the 
CCFULL calculation considering the two-phonon coupling. 
It shows a part loss of the component below 75 MeV. This is 
due to the coupling to the Q >0 neutron transfers, correspond- 
ing to the further fusion enhancement at sub-barrier energies 
compared with the calculation. 
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FIG. 6: Barrier distributions for ^"S+'^Zr from fusion (solid circles) 
and quasi-elastic scattering (hollow circles). The line is the CCFULL 
calculation with two-phonon coupling. 
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FIG. 7: For ^^S+'*Zr, the symbols are the same as in Fig. 6. 



V. SUMMARY 

The fusion excitation functions for ■^^SH-''"'^^Zr were mea- 
sured with a high precision near and below the Coulomb bar- 
rier. The sub-barrier cross sections for ^^S-H^^Zr are much 
larger compared with ^^S-n^^Zr The data have been analyzed 
with both the CCFULL code and the NTFus code based on 
a semi-classical coupled-channels model, which includes the 
sequential neutron transfers for ^^S+^*Zr as earlier proposed 
by Zagrebaev. Good agreement between experimental data 
and the calculation is achieved for ^^SH-''"Zr by including the 
couplings to the low-lying quadruple and octupole vibrations 
in ^-S and '"'Zr. Whereas the NTFus calculations can repro- 
duce the data by including four sequential neutron transfer 
channels as well as the low-lying quadrupole and octupole 
vibrations in ^^S and '"'Zr. The comparison with previous 
data on '**'Ca-i-^^Zr shows that the excitation functions of ■'^S, 
'"'Ca+^^Zr are very similar in the whole energy range. Both 
systems have similar collective states and positive Q-valne 
neutron transfer channels. This comparison strongly supports 



the previous suggestion fls] that positive Q-value neutron 
transfer channels enhance sub-barrier fusion cross sections, 
particularly at very low energies. Also the fusion and quasi- 
elastic barrier distributions of the ■'^S+^^'^^Zr systems are ba- 
sically consistent in both cases. The two barrier distributions 
for ^^S-H^^Zr cannot be reproduced by CCFULL code below 
75 MeV. The fact shows again the effect of the Q >0 neu- 
tron transfers on the sub-barrier fusion process. In addition to 
the fusion excitation function, the neutron transfer cross sec- 
tion measurement for this system should provide useful infor- 
mation on the coupling strength of neutron transfer channels, 
which will allow us to reach a much deeper understanding of 
the role of neutron transfer mechanisms, sequential or simul- 
taneous, in the fusion process. 
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